Mechanical Response of a Small Swimmer Driven by Conformational Transitions 
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A conformation space kinetic model is constructed to drive the deformation cycle of a three-sphere 
swimmer to achieve propulsion at low Reynolds number. We analyze the effect of an external load 
on the performance of this kinetic swimmer, and show that it depends sensitively on where the force 
is exerted, so that there is no general force-velocity relation. We discuss how the conformational 
cycle of such swimmers should be designed to increase their performance in resisting forces applied 
at specific points. 
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Active transport is a most fascinating aspect of the 
busy life in the cell [l|. In the nano-scale world where 
thermal agitations are wild, miniature machines called 
molecular motors convert chemical energy — from hydrol- 
ysis of ATP molecules — directly into useful mechanical 
work, in the form of carrying cargo or sliding actin fila- 
ments along one another. While it is difficult to imagine 
fabricating such sophisticated machines in the lab, one 
may naturally wonder if it is possible to design simpler 
machines with similar functionalities [2|. In this flavor, 
an interesting target is an autonomous small scale swim- 
mer 0, 1J] , which could later on be steered by coupling to 
a guiding network or system. 

Swimmers at small scale (low Reynolds number) have 
to undergo non-reciprocal deformations to break the 
time-reversal symmetry and achieve propulsion [5| . This 
imposes significant constraints when one wants to design 
a swimmer with only a few degrees of freedom and strike 
a balance between simplicity and functionality (a]- Re- 
cently, there has been an increased interest in such de- 
signs [71, y and an interesting example of such robotic 
micro-swimmers has been realized experimentally using 
magnetic colloids attached by DNA-linkers [3|. 

Here we combine features of simple low Reynolds num- 
ber hydrodynamic swimmers and elements characteristic 
of models for chemical molecular motors. We focus on 
a recently introduced three-sphere swimmer [7| with the 
minimal two degrees of freedom. Instead of assuming a 
prescribed sequence of deformations, we consider these 
deformations to occur stochastically, as conformational 
transitions between elongated and shortened states for 
each of the two degrees of freedom. This gives us a swim- 
mer with a velocity that depends on the transition rates 
between these states, which in practice could come about 
via mechanochemical transitions, i.e. due to chemical 
reactions that are coupled with such mechanical defor- 
mations [1[. We check that a net velocity requires that 
detailed balance in the transition rates is broken. Us- 
ing this simple kinetic model, we study the effect on the 
swimming velocity of a resisting external force or load: 
the load clearly drags the swimmer backwards, but also 



puts elements in compression or extension, thereby mod- 
ifying the transition rates between extended and short- 
ened states. As a consequence, we find that the perfor- 
mance of the motor strongly depends on where the force 
is exerted, in contrast to the usual perception that the 
performance of a swimmer - or a motor - can be summa- 
rized in a unique force-velocity relation. Interestingly, 
the motor performance can in some special cases be in- 
creased upon application of the external load, provided 
it is applied at the right location. More generally, we dis- 
cuss efficient strategies for optimizing the performance of 
this swimmer. 

We start with the swimmer model introduced in [7|, 
namely three spheres connected by two linkers of neg- 
ligible hydrodynamic effect, that cycle in time between 
extended and short states. For simplicity, we take here 
all sphere radii to be equal to a. We further assume 
that the lengths of the two arms are Li{t) — ii + ui{t) 
and L2{t) = £2 + "2(0 with the Ui's being small per- 
turbations about the average lengths. For prescribed 
arm deformations, writing a force balance on each sphere 
leads an instantaneous net displacement velocity of the 
swimmer, that can here be written as a series expansion 
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FIG. 1: Schematics of the configuration space of the three- 
sphere swimmer. To maintain a net swimming to the right, 
the deformation moves need to make more clockwise full cy- 
cles than counterclockwise ones. 
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FIG. 2: Force balance on the spheres and the hnkers for a 
swimmer moving to the right, when the external load F (act- 
ing to the right) is attached to the (a) head, (b) tail, or (c) 
middle of the swimmer. In each case, it is identified whether 
each linker is under compression (C) or tension (T), with the 
value of the force given (underneath). 



vit) = AiUi + BijUiUj+CijkUiUjUk + - ■ ■ , where the coef- 
ficients Ai, Bij, Cijk, etc. are purely geometrical prefac- 
tors (i.e. involving only the length scales a and £i). After 
many cycles, this process gives a vanishing contribution 
from the linear terms ui and U2 and from the symmetric 
combination uiU2 + U2W1 = d{uiU2)/dt. Thus to leading 
order the average swimming velocity is 
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where dA is the area element in the (^1,^2) space, and 
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J9|. In other words, to the 
leading order the swimming velocity is proportional to 
the area enclosed by the orbit of the cyclic motion in the 
configuration space of the deformations. 

We now focus on a situation where the two arms can 
be in two states with deformations of either m^ = or 
Ui = 5i, and transit from one to the other in an almost 
instantaneous fashion. This means that the configuration 
space of the swimmer will be made of only four distinct 
states as shown in Fig. [l] which correspond to different 
values of the pair {ui,U2), namely: state A for {5i,82), 
state B for ((5i,0), state C for (0,0), and state D for 
(0,(52). We then assign transition rates to the system, 
corresponding to the average rate of opening and closing 
of the arms. For example, the transition rate from state 
A to state B is denoted as kg a, and similar notations 
are used for the 8 rates describing forward and reverse 
transitions along the cycle 
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For simplicity, and at the cost of motor efficiency, we as- 
sume that the transitions occur quite rapidly and seldom, 
so that they never "overlap." 

We can now calculate the swimming velocity as a func- 
tion the transition rates. At steady state, the average 
swimming velocity of the object is given by the probabil- 
ity current J along the A— >B— >C^D— s-A cycle times the 
net displacement while performing the cycle. This dis- 
tance Ax is simply K5i52, which yields V = K6i52J. The 
probability current J is a function of the transition rates, 
which can be obtained from straightforward algebra: 



J = 
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From the above equation it is clear that if detailed 
balance holds, then J is zero as the numerator van- 
ishes. Using the average steady state current, we can 
deduce the average period of one full cyclic motion along 
A^B^C^D^A as T = J-i. In general a 1 < — > 2 
asymmetry in the system, together with breaking the de- 
tailed balance at least for one of the transitions, will 
lead to net motion. For the particular limit where 
the forward rates are all much higher than the corre- 
sponding backward ones (fcsA ^ kAB, etc.), we find 



T = k in + kjj'jj + fcpjj -j- kg^, which simply means that 
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the period for a full cycle is the sum of the time inter- 
vals needed to complete each leg of the cycle. As an- 
other example, we can assume that all of the equilibrium 



kfSaS are equal to 1 (for the sake of illustration), and 
that by external action only one of them is modified e.g. 
kBA = 1 + e. In this case, it is easy to show that Eq. 
([3]) yields J = e/(16 -f 6e), which leads to a velocity pro- 
portional to the perturbation if the latter is small and 
independent of it if the perturbation is very large, as the 
cycling is then limited by the other three unperturbed 
transitions. In general, it is easy to see that the slowest 
leg of the reaction controls the average rate of full cyclic 
motion. 

It is interesting to study the effect of an external load 
on the velocity of the system and the performance of the 
motor. When the swimmer is subject to external forces, 
because of carrying a cargo for example, there are two 
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FIG. 3: Motor performance function J{F)/J{0) versus the external force, when only one of the rates is enhanced to 1 + e while 
the rest are kept fixed at 1. The plots correspond to the enhanced rate being (a) koBA, (b) fcocs, (c) koDC, and (d) koAD, 
and e = 50. In each case, the dashed (red) line corresponds to the load being attached to the head, the dotted (blue) line 
corresponds to attachment to the tail, and the solid (green) line is for the middle attachment. 



types of mechanical responses in the system. First, the 
external forces enter the hydrodynamic force balance on 
each sphere, and this will introduce a Stokes drag on the 
sphere as a whole, which is a linearly decaying contribu- 
tion to the net swimming velocity as a function of force. 
Second, the transition rates that control the kinetics of 
the deformations for the two arms of the swimmer arc 
affected by the external forces as they will have to do 
mechanical work against them to induce the deforma- 
tions. Depending on where the force is applied, different 
legs of the kinetic cycle could be affected, and this could 
lead to a complex mechanical response with the perfor- 
mance of the motor depending on the location of the load. 
The force-dependent kinetic rates will yield a net current 
J{F), which combines with the Stokes response to give 
the swimming velocity as 
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where Vq is the swimming velocity at zero force, rj is 
the viscosity of the solvent, and a/? is a rcnormalizcd 
hydrodynamic radius [10|. 

The transition rates are modified in the presence of ex- 
ternal forces, because the mechanical energy enters the 
balance of probability of the different states and transi- 
tions among them. If there is a transition from a — > /? 
that corresponds to an extension by a factor of 6, then 
under a positive tension / the rate of a — > /3 transi- 
tions is increased by a factor of exjp{fx/k^T) while the 
rate of /3 ^ a transition is decreased by a factor of 
exp(— /cc'/fceT) where typically x = 66 \s the distance 
between state-a and the energy barrier and x' = (1 — 9)5 
is the distance between the energy barrier and state-/3 
{9 between and 1). Thus the ratio between the two 
transitions rates [i.e. the a ^ (3 rate divided by the 
/? — > a rate) changes under a tension / by a factor 
of exp(/d/fcBT) as required by the Boltzmann formula 
(equilibrium populations between state-a and state-/? un- 
der /). 

In our system, the value of the force under which each 



arm should close or open depends on where the load is 
applied. Figure [2] shows the break down of the mechani- 
cal force balance on each sphere, and the corresponding 
forces endured by each linker, for the three different po- 
sitions of the load. When the resisting force is attached 
to the head of the swimmer, both linkers are under com- 
prcssional forces, and the compression force on the right 
arm — nearer to the load — is larger than that of the left 
arm by a factor of two. Attaching the load at the tail cre- 
ates a similar pattern of tensional forces. If the force acts 
on the middle sphere, the left arm is under compression 
and the right arm is under tension. 

Using the above definition, the transition rates from a 
conformation state a to another state /3 can be written 
as 



fc/3a = kopa exp 



1 f/3aSi 

2 /cbT 



(5) 



where ff3a is the force endured by the linker i that un- 
dergoes a deformation during the a — > /3 transition, and 
9 = 1/2 is assumed for simplicity. The sign of fpa is 
determined by whether the transition (deformation) is 
helped {+) or opposed (— ) by the force acting on the 
linker. The values of fi3a are given in Table U for the 
forward reaction rates for the different locations of the 
load. Note that by definition ffSa = —faf3, which can be 
readily used to calculate the reverse rates. 

The force-dependent rates [from Eq. (O and Table [I] 



TABLE I: The algebraic force fpc, which should be used in 
Eq. ([5} to calculate the forward rates. The values for the 
corresponding reverse rates can be obtained via fp^ = —/a/3. 
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can be used in Eq. Q to calculate the current J{F), 
which determines the swimming velocity under the effect 
of an external load F. From Eq. ^, it appears that the 
normalized current J{F)/J{0) is a quantitative measure 
of how the ability of the motor to generate propulsion is 
affected by the presence of the load. In Fig. [31 this "mo- 
tor performance function" is plotted against the external 
force for the particular example discussed above (in which 
only one of the forward rates is enhanced to 1 + e while 
the rest of the rates are set to unity), and 5i ~ 82 = S 
is assumed for simplicity. Figure [3^ corresponds to when 
the A ^ B (contraction of the left arm) transition rate 
is enhanced, and it shows that attaching the load to the 
head or the middle for both of which the left arm is under 
a compression of ^F quickly decreases the performance 
of motor. On the other hand, attaching the load to the 
tail of the swimmer, which puts the left arm under a ten- 
sion of |F, actually helps the motor initially for forces 
of up to ik^T /5 or so, before eventually hampering the 
performance at large forces. One notes that the force 
across the left arm actually helps the A ^ B transition 
when the load is at the head or the middle, and opposes 
it when it is at the tail. It thus seems that the perfor- 
mance of the motor is best when the rate is enhanced for 
the deformation which is most hampered by the exter- 
nal load. In other words, the best strategy seems to be 
to try and make the performances of the different legs 
of the reaction cycle as uniform as possible, as the total 
velocity is controlled by the weakest performance in the 
cycle. The same pattern can be seen in Figs. [3|d-[311. An- 
other interesting feature that can be seen is that when 
the load is at the middle and the condition is right for 
improved performance (see above) the system seems to 
endure comparatively much stronger forces: in Figs. [8)3 
and [31; one can sec that the performance is significant 
for loads of up to about l2kBT/S. This is presumably 
because attaching the load to the middle creates a more 
balanced distribution of the forces in the linkers (still of 
opposite nature but of equal magnitudes; see Fig. [2land 
Table IB). 

Even when the performance of the motor is increased 
by the opposing force, one still has a decreasing trend 
for the swimming velocity because of the Stokes drag 
term in Eq. ([4]). Using a linear approximation for 
J{F) ~ Jo(l + cFS/kBT) at small forces (where c is a 
positive constant of order unity), one can write Eq. (|4]) 



as V{F) = Vo 
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, which implies 



187r?7a/^Vo 

that for forces much smaller than the thermal activation 
force kBT/S the increased motor performance can lead 
to increased swimming velocity if the viscous drag on 
the swimmer is larger than the thermal activation force. 
While this could be extremely difficult to achieve as it 
requires unrealistically high swimming velocities, it is an 



interesting fundamental possibility that increased swim- 
ming velocity can be achieved upon exerting opposing 
forces. 

In conclusion, we have proposed and studied a simple 
model of a low Reynolds number swimmer driven by a 
kinetic engine. The main result is that the ability of this 
swimmer to carry a load or to resist an opposing force 
depends on where the load or the force is applied. This is 
not linked to the stochastic nature of the present motor, 
but also holds for a motor driven by a prescribed sequence 
of internal stresses (to which the applied stresses add 
up). Altogether, this shows that the description of such 
machines can go beyond a simple force-velocity relation, 
more complex and maybe richer in functionality. 
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■J-), which should be expanded in 2nd order of 
the deformations to keep the same order of terms in the 
perturbation theory. 



